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1.  INTRODUCTION 


A.  BACKGROUND 

Dynamics  play  a  major  role  in  the  design  and  development  of  the  modern-day 
helicopter.  Helicopter  vibratory  characteristics  have  been  the  decisive  factor  in  the 
selection  of  winning  designs  in  helicopter  program  competitions,  as  well  as  a  major  cause 
of  cancellation  of  other  programs.  Low  vibratory  levels  increase  fatigue  life  of 
components,  improve  passenger  comfort,  and  decrease  crew  fatigue.  Therefore,  it  is  of 
utmost  importance  to  be  able  to  analyze  and  examine  the  effect  of  these  vibrations  in  the 
preliminary  stages  of  helicopter  design. 

In  helicopters,  vibrations  fall  into  three  categories:  vibrations  due  to  rotor  excitation, 
which  occur  at  frequencies  that  are  integral  multiples  of  the  rotor's  rotational  speed; 
vibrations  due  to  random  aerodynamic  excitation,  where  the  frequencies  observed  are  the 
natural  frequencies  of  the  structure  being  excited,  and  self-excited  vibrations  such  as 
flutter  and  ground  resonance.  In  this  paper  the  primary  source  of  helicopter  vibrations  will 
be  considered,  that  is,  the  rotor  system. 

The  alternating  forces  that  excite  the  rotor  blade  are  almost  entirely  due  to  the 
periodic  variations  in  blade  airloads  encountered  in  flight.  These  time-dependent 
aerodynamic  forces  occur  due  to  the  variations  in  velocity  and  angle  of  attack 
encountered  by  the  rotating  blade.  In  low-speed  flight  blade-vortex  interaction  plays  a 
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significant  role  in  blade  angle  of  attack  variation  due  to  the  interaction  of  a  rotor  blade 
with  the  trailing  tip  vortex  shed  by  the  preceding  blade.  In  steady-state  level  fight  these 
loads  can  be  considered  harmonic  since  they  occur  at  multiples  of  the  rotor's  rotational 
speed,  i.e.  1/rev.,  2/rev.,  3/rev.,  etc.  The  blade's  ensuing  response  to  these  loads  is 
transferred  as  vibratory  shears  and  moments  to  the  rotor  head  and  filter  through  to  the 
fuselage  Therefore,  the  ability  to  analyze  the  stresses  in  the  helicopter's  rotor  blade  will 
lead  to  an  understanding  of  how  to  alleviate  these  stresses  To  that  end,  this  thesis 
presents  a  computational  tool  to  examine  the  flatw  ise  dynamic  loads  of  a  rotor  blade  with 
vaiying  material  properties 
B.  JANRAD 

The  Joint  Army/Nav'y  Rotorcraft  Analysis  and  Design  (JANRAD)  computer 
program  was  developed  to  meet  the  speciflc  needs  of  the  Naval  Postgraduate  School 
(NPS)  for  preliminary  helicopter  design.  JANRAD  was  written  as  an  interactive,  user 
friendly  program,  capable  of  accurately  and  quickly  solving  helicopter  design  problems  at 
the  preliminary  design  level.  JANRAD  consists  of  three  major  subroutines.  The  f  rst 
routine,  JANRAD  Performance,  calculates  the  trim  solution  and  various  performance 
parameters  used  in  helicopter  design.  This  program  is  described  in  detail  in  Ref  1 .  The 
second  component  is  JANRAD  Stability.  This  routine  calculates  stability  derivatives,  and 
describes  the  open  loop  control  characteristics  using  a  linear  state  space  model.  The 
details  of  this  program  are  provided  in  Ref  2. 
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The  third  component  is  JANRAD  Dynamics,  written  by  the  author,  which  adopts 
the  same  user  friendly  menu  driven  format  to  assist  in  determining  preliminary  rotor  blade 
design  parameters.  The  program  uses  a  modified  Myklestad-Prohl  Method  [Ref  3]  to 
determine  the  uncoupled,  flatwise,  forced  response  of  the  helicopter  rotor  blade  to 
aerodynamic  loads  and  dynamic  forces.  These  forces  are  calculated  using  the  same 
helicopter  design  parameters  and  data  output  from  the  JANRAD  Performance  routine, 
making  JANRAD  a  truly  integrated  design  tool. 

JANRAD  was  originally  written  for  version  3.5  of  386-MATLAB*  but  since  has 
been  updated  to  MATLAB*"  version  4.0  for  WINDOWS’^.  MATLAB  is  a  high 
performance  interactive  software  package  for  scientific  and  engineering  numeric 
computations.  This  program  was  chosen  due  to  its  wide  accessibility  at  the  NPS  and  for 
its  ease  of  use.  This  frees  the  engineering  student  from  the  burden  of  learning  "low  level" 
programming  languages  and  allows  the  user  to  concentrate  on  the  design  process  itself 
Athough  there  are  other  helicopter  design  programs  such  as  HESCOMP,  CAMRAD  or 
RACAP;  these  programs  are  often  proprietary,  expensive  to  acquire,  and  cumbersome  to 
use.  Whereas  the  general  simplicity  of  MATLAB'^  allows  the  NPS  student  to  modify 
JANRAD  to  meet  his/her  specific  needs. 
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11.  THEORY 


There  are  several  methods  for  determining  the  forced  response  of  dynamic  systems. 
These  methods  may  be  classified  into  two  categories:  exact  methods  and  approximate 
methods.  In  exact  methods  an  explicit  solution  of  the  governing  differential  equation  is 
obtained  in  closed  form.  In  general,  explicit  solutions  are  possible  only  for  linear 
differential  equations  with  constant  coefficients.  These  types  of  equations  are  associated 
with  small  amplitude  vibrations  of  structures  having  uniformly  distributed  mass  and 
stiffness  properties.  When  these  properties  are  not  uniformly  distributed,  the  coefficients 
are  variable  and  exact  solutions  are  not  always  possible.  These  problems  in  transverse 
vibrations  of  beams  can  be  solved  by  an  approximate  numerical  procedure  which  was 
developed  independently  by  Prohl  [Ref  4]  and  Myklestad  [Ref  5].  As  a  historical,  as  well 
as  a  pedagogical  precursor  to  the  Myklestad-Prohl  method,  the  Holzer  method  [Ref  6, 
Ref  7]  demonstrates  a  simpler  lumped  mass  system  that  models  torsional  dynamics. 

A.  HOLZER'S  METHOD  FOR  TORSIONAL  SYSTEMS 

The  Holzer  method  is  an  approximate  solution  to  the  second  order  differential 
equation  of  a  uniform  slender  shaft: 


where 


dx^ 


+ 


co^/ 

GJ 


0  =  0 


0  =  maximum  amplitude  of  torsional  vibration 


(1) 
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co=  frequency  of  a  natural  mode 
I  =  mass  moment  of  inertia  per  unit  length 
GJ  =  torsional  stiffness 

The  Holzer  method  approximates  the  differential  equation  of  motion  by  two 
difference  equations  which  replace  the  distributed  structure  with  one  in  which  the  masses 
and  connecting  elastic  elements  are  discretized  or  "lumped."  These  numerical  models  are 
often  referred  to  as  discrete  or  lumped  mass  sy’stems.  Figure  1  shows  a  part  of  the 
non-uniform  shaft,  lumped  mass  system  used  in  deriving  the  Holzer  difference  equations. 


Figure  1.  Ideal  Lumped  Mass  System  For  A  Non-Uniform  Shaft  From  Ref.  7 
The  equation  of  moment  equilibrium  for  the  rth  mass  is: 


Mfi+i  ~  Mfi  +  =  0 


(2) 
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where 


,  =  torque  applied  to  the  mass  /  +  1  by  shaft  section  i 
=  torque  applied  to  the  mass  /'  by  shaft  section  /  -  1 
1  =  moment  of  inertia  of  mass  / 

i 

0,  =  angular  displacement  of  mass  / 

03  =  angular  frequency  of  vibration 

A  second  equation  is  derived  by  expressing  the  twist  of  the  shaft  section  i  in 
terms  of  the  torque  on  that  section. 

e„,  -  e,  =  ^  (3) 

ki 

where 

k,  =  stiffness  of  shaft  section  i 

Equations  2  and  3  can  be  written  in  the  following  form  for  purposes  of  an  iterative  analysis 


Mti+i  —  ~  O)^//0j 


(4) 


Ji+i 


0,  + 


Mfi+i 

ki 


(5) 


Equations  4  and  5  may  also  be  derived  from  the  differential  equation  for  a 
uniform  slender  shaft.  Equation  1.  The  parametric  form  of  Equation  1  may  be  written  as 


^=-a)V(x)0  (6) 

ax 
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and 


(7) 


cJQ  Mr 
dx  G./(.v) 


If  the  first  derivatives  are  approximated  by  the  difference  forms  as 


d A  / y  ^  A  f  i  i+\  A  / 1  I 

dx  ~  Ax. 


d%  ^  e,.i  -  e, 

dx  A.v, 


and  bv  usiim  the  definitions 


/,  =  /(.v)Av, 


GJix) 


we  may  once  again  obtain  the  Holzer  difference  Equations  4  and  5. 

B.  MODIFIED  MYKLESTAD-PROHL  METHOD 

An  extension  of  the  Holzer  method  for  torsional  systems  is  the  modified 
Myklestad-Prohl  method  for  lumped  mass  beam  systems  [Ref  8],  While  the  Holzer 
method  solves  a  second  order  linear  equation  with  two  difference  equations,  the 
Myklestad-Prohl  method  solves  the  fourth  order  differential  equation  for  a  uniform  beam: 


r,jd^Z  (fz 

tl — 7  =—m  — - 
dx^*  di^ 
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with  four  difference  equations,  where 

Z  =  deflection,  expressed  in  terms  of  vertical  displacement 
dx  =  spanwise  differential  length 
m  =  beam  mass  per  unit  length 
El  =  beam  stiffness  or  bending  modulus 

The  modified  Myklestad  method  uses  these  four  difference  equations  to 
progressively  compute  the  shear,  moment,  slope,  and  deflection  from  one  section  to  the 
next,  taking  the  centrifugal  and  aerodynamic  forces  into  account.  Figure  2  illustrates  an 
element  of  the  blade  lumped  mass  model  with  the  uncoupled  flatwise  dynamic  forces  and 
moments  acting  on  it. 


Figure  2.  Blade  Element  Lumped  Mass  Model  For  Uncoupled  FlaUvise  Vibration  From  Ref.  3 
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By  taking  the  blade  element  model  into  consideration  it  is  possible  to  write 

equations  for  the  tension,  shear,  and  moment  at  //  entirely  in  terms  of  quantities  at  //  1. 

These  can  then  be  substituted  into  the  geometric  equations  for  the  slope,  0  and  the  \'ertical 
deflection,  Z.  Writing  the  equations  of  equilibrium  for  the  //th  element  of  the  uncoupled 

flatwise  system  we  ha\’e: 

Centrifugal  Force: 

in  —  N+\  n  ^  ^  n  (  k"*! 

where; 

T  =  centrifugal  tension  force  acting  on  the  //th  blade  element 

ny^  =  mass  of  //th  blade  element 
Q  =  rotor  angular  velocity 

q,  =  distance  of  the  //th  blade  element  from  the  axis  of  rotation 

Shear: 

Sn  —  iS  fi  +  l  +  HI  t)  (Si  7.  n  j(  n  (Si  7  n  +  /;  +  j.f  n  (^‘^) 

where; 

S,  =  flatwise  shear  force  acting  on  the  //th  blade  element 

m,,  =  mass  of  //th  blade  element 
(0  =  applied  excitation  frequency 

=  deflection,  expressed  in  terms  of  vertical  displacement  of  the  //th  blade 
element 

kn  ^  aerodynamic  thrust  acting  on  the  //th  blade  element  for  a  particular 
frequency  component 

C,,  =  flatwise  aerodynamic  damping  for  the  //th  blade  element 
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The  aerodynamic  damping  term  can  be  expressed  by: 


C„  =  -7-^  (c/7or<y) „(/„,„+])  (p/2)  (Q/-„) 
da 

where; 

=  lift  coefficient  slope  for  two-dimensional  blade  section 

da 

/  ,  =  length  of  »th  blade  element 

p  =  mass  density  of  air 
Q.  -  rotor  angular  velocity 

r  =  distance  of  the  mh  blade  element  from  the  axis  of  rotation 

n 

Additionally,  from  Figure  2,  we  have  the  equations  for  the  moment: 


where: 


M n  —  M f]+\  T  S n+\  I n.n+\  ^n+I  (^n+]  Z n) 


=  flatwise  bending  moment  at  rrth  blade  station 
=  flatwise  shear  force  acting  on  the  wth  blade  element 
n+]  -  length  of  «th  blade  element 

=  centrifugal  tension  force  acting  on  the  I  blade  element 
=  deflection,  expressed  in  terms  of  vertical  displacement  of  the  »th  blade 
element 


slope: 


6  n  -  0  n+l  \  +  T 
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where; 

0,;  =  flatwise  slope  of  the  elastic  axis  at  the  //th  blade  station 
=  centrifugal  tension  force  acting  on  the  //th  blade  element 
/  ,  =  lenizth  of  //th  blade  element 

El  =  beam  bending  modulus  at  the  //th  blade  station 
M  =  flatwise  bending  moment  at  //th  blade  station 
S  ,  =  flatwise  shear  force  acting  on  the  //  1  blade  element 

n- I  ^ 


and  deflection; 

Z„  =  Z„+i 

where; 


6  n  /  n. 


H+l 


t  n+l  6«+l  ^  n.n+] 

JE!  iW 


n+l 


/2 

'  n.  n+l 


s 


«+l  '  n.n+1 


3  El 


(18) 


Z,,  =  deflection,  expressed  in  terms  of  vertical  displacement  of  the  //th  blade 
element 

0  =  flatwise  slope  of  the  elastic  axis  at  the  //th  blade  station 
/„  =  length  of  //th  blade  element 

=  centrifugal  tension  force  acting  on  the  //th  blade  element 
El  =  beam  bending  modulus  at  the  //th  blade  station 
M,,,  I  =  flatwise  bending  moment  at  //  '^  I  blade  station 

=  flatwise  shear  force  acting  on  the  1  blade  element 
The  equations.  Equations  14,16,17,  and  18,  may  be  rewritten  so  that  all  the 
terms  on  the  right  side  of  the  equation  are  in  terms  of  //+ 1  variables.  By  then  expressing 
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these  equations  in  matrix  form  and  integrating  the  blade  elements  from  tip  to  root ,  the 
shear,  moment,  slope,  and  deflection  at  the  blade  root  may  be  expressed  in  terms  of 
corresponding  unknown  tip  values.  With  the  use  of  the  appropriate  root  boundary 
conditions  the  unknown  tip  values  may  be  determined.  Once  this  is  done  the  shear, 
moment,  slope,  and  deflection  can  be  determined  along  the  length  of  the  rotor  blade.  The 
details  of  this  procedure  will  be  discussed  in  the  Blade. m  program  details  section. 
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Ill.  M-FILE  DETAILS 


A.  DYNAM.M 

The  first  module  of  the  blade  dynamic  analysis  portion  of  JANRAD  is  Dynam.m. 

The  primary  function  of  this  subroutine  is  to  enter  and  store  rotor  blade  data  required  for 
analysis.  Since  this  portion  of  JANRAD  uses  variables  which  are  cumbersome  to  keep 
re-entering,  it  also  follows  a  similar  menu  driven  format  that  allows  the  user  to  save 
variables  for  later  analysis.  Dynom.ni  uses  output  data  from  the  JANRAD  Rotor 
Performance  program  so  it  must  first  determine  whether  it  has  been  run.  If  it  determines 
that  the  Rotor  Performance  has  not  been  run,  the  user  is  exited  out  of  the  program  to 
enter  the  JANRAD  main  menu  again.  The  program  then  asks  the  user  to  either  edit  or 
create  a  data  file.  If  editing  is  chosen,  the  program  then  asks  the  user  for  the  file  name.  If 
the  file  is  not  found,  the  user  is  instructed  to  try  again  or  exit  the  program.  When  the  file 
is  found,  an  edit  menu  is  presented  where  the  user  may  change  the  rotor  blade  parameters. 

Whether  the  user  chooses  to  create  a  new  data  file  or  to  edit  a  previous  one,  the 
program  prompts  the  user  to  enter  the  rotor  blade  material  properties  and  weight 
distribution.  The  most  frequently  applied  boundary  conditions  for  the  root  end  of  the 
rotor  blade  are  either  articulated  or  hingeless.  For  the  articulated  case,  the  flatwise 
moment  is  zero,  whereas  in  the  hingeless  case  the  root  slope  is  equal  to  zero.  In  both 
cases  the  blade  root  vertical  displacement  is  zero. 
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The  material  properties  may  he  entered  in  two  ways,  either  as  a  combination  of  the 
modulus  of  elasticity,  E  (Ib/in"),  and  the  rotor  blade  flapwise  moment  of  inertia,  I,,  (irE),  or 
in  combination  as  the  rotor  blade  stiffness  coefficient,  El  (lb-in").  These  variables  along 
with  the  weight  distribution  are  entered  as  vectors  whose  elements  represent  the  elemental 
properties  along  the  blade.  The  user  is  instructed  to  enter  the  rotor  blade  properties  from 
blade  tip  to  blade  root  and  is  reminded  to  enter  the  same  number  of  blade  elements  as  was 
previously  chosen  in  the  main  JANRAD  menu.  When  entering  the  weight  distribution,  the 
weight  is  given  as  previously  entered  for  the  Rotor  Performance  portion  of  JANRAD. 

The  entered  data  vectors  are  then  checked  for  correct  length  and  stored.  The  data  is 
saved  using  a  file  name  in  the  same  manner  as  used  in  the  main  JANRAD  program.  These 
specifics  are  outlined  in  Ref  1 . 

B.  BLADE.M 

Blade,  m  uses  the  Mylestad-Prohl  equations  for  rotating  beams  in  complex  form  to 
calculate  the  flatwise  forced  response  of  a  rotor  blade  of  var>'ing  material  properties.  This 
process  is  accomplished  in  four  steps: 

1 .  Calculate  the  aerodynamic  damping,  the  centrifugal  force,  and  the  steady  and 
first  ten  harmonic  thrust  differentials. 

2.  Calculate  the  transfer  matrices  for  the  harmonic  blade  response  along  the 
length  of  the  blade. 

3.  Solve  for  the  rotor  blade  tip  unknowns  (Z,,  6?  ). 

4.  Repeat  steps  3  and  4  utilizing  the  respective  thrust  harmonic  and  excitation 
frequency. 
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Before  discussing  the  four  step  process  used  in  the  program,  we  must  return  to  the 
Myklestad-Prohl  equations  and  present  them  in  the  form  they  are  used  in  Blade. m. 

The  modified  Myklestad  equations,  which  were  introduced  earlier  as  Equations 
14,  16  -  18,  may  be  rewritten  so  that  all  the  terms  are  on  the  right  side  of  the 
equations,  as  follows: 


S„  =  d\\Sn  +  \  +  d]2M„+]  +  d  13  Qn+\  +  d  1 4  -Z„+l 


M„  =  d2\S„+\  +  d  22  M  n+\  +  d  23^  n+\  +  d  24  ^ 


0,,  =d3\Sn4\  +d32M„4.\  +6^33  6„+l  +d34Zn4\ 


Zn  =  dA\S  n+X  +i^42-^n+l  d  43%  „  +  \  +  d^n  Zn+X 


where  d  between  station  n  and  station  w+ 1  is  given  by: 


d\x  =  1  +(WnCjC)^  -yC„(o)- 


f  ('23') 


f  T  B 

d  \3  —  —(r)1„03  —jCnOj)  ()E1  '«,«+! 


di4  =  -jCrtOi)  (24) 


,  T„  +  \  l„,n+X  ,  ,  _  .4 

"21  =  - - 1"  'n,«+l  -  “"43 

6L1 


,  Tn+ll  «.«+!  .  1  J  /oc\ 

d22  =  - YeI -  +  1  =  ^33  (25) 


d23  —  — 


rp2  p 

‘  n  +  X  '  «,«+! 


"I"  T n+\  I n,n  +  l 


d24  —  0  —  d3 
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1  ’  nji-T  1  I 

‘^31  -  -  2^r/  “  “ 


c/  3  2  ■“ 


c/41  — 


'  /..//+ 1 


These  equations  nia\’  also  be  written  in  the  matrix  form 


where 


!Cl  „  +  i  + 


tCK,  = 


1 1 

dM 

d\y 

d  1  4 

2  1 

dll 

dll 

di4 

1 

dii 

dll 

di4 

4  1 

d  Ai 

d  41 

c/44 

Matrix  Equation  29  relates  the  response  of  the  // 1  1  element  to  the  adjacent  inboard  //th 
element  through  the  transfer  matrix  It  should  be  noted  that  [D]^,,  is  dependent 

only  on  known  or  pre-calculated  parameters.  The  elements  in  Equation  32  represent  the 
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shear,  moment,  slope,  and  deflection  which  are  solely  due  to  aerodynamic  loads.  Now 
that  the  Myklestad-Prohl  equations  have  been  presented  in  the  matrix  form  used  in  the 
program,  we  will  proceed  to  discuss  the  first  procedure  of  Blade.m. 

The  flatwise  aerodynamic  damping  term,  C„,  for  the  mh  term  is  given  by: 

^  (chord) „  )  (p/2)  (Qr„)  (33) 

da. 

where 

=  lift  coefficient  slope  for  two-dimensional  blade  section 
da 

I  ,  =  length  of  //th  blade  element 
p  =  mass  density  of  air 
Q  =  rotor  angular  velocity 

r  =  distance  of  the  «th  blade  element  from  the  axis  of  rotation 

n 

The  lift  coefficient  slope,  — as  well  as  the  mass  density  of  air,  p,  and  the 

da 

rotor  angular  velocity,  Q,  have  been  previously  entered  as  required  variables  for  JANRAD 
rotor  performance.  The  elemental  chord  length  has  been  updated  to  account  for  taper 

ratio  and  is  also  taken  from  JANRAD. 

In  calculating  the  steady  and  first  ten  harmonic  thrust  differentials  the  elemental 
thrust  forces  along  the  blade  are  obtained  from  the  rotor  performance  analysis  portion  of 
JANRAD.  Using  a  numerical  Fourier  analysis,  these  thrust  forces  are  resolved  into  a 
steady  and  ten  harmonic  components.  These  forces  may  be  expressed  in  complex  form  as; 

10 

Fx,;(r„)  =  Fo(r„)  +  X  %eal(fi{r„)-jFi(r„))e^"'^  (34) 

i=\ 
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where  the  steady  aerodynamic  thaist  is  given  by: 


/■ 


and  the  harmonic  thaist  components  are: 


/*'/(/■-;)  =  77-  Z  r/7'v(/„)  sin(/V|/)  /  =  1,2  . .10 


f,{rn)  =  jf;  Z  dl\/{r„)  cos {i\f)  i  =  1,2... 10 

j=  I 


where: 


n  y  =  number  of  azimuth  sectors 

d'[\,(r„)  =  elemental  thaist  at  blade  azimuth  position  y,  and  blade  station  r^ 
y  =  azimuth  angle  position 
/  "  /th  harmonic  component 

The  second  step  in  Blade. ni  involves  progressively  integrating  Equation  29  from 
tip  to  root.  The  general  form  of  Equation  29  may  be  obtained  prior  to  substituting  the 
boundary  conditions  for  either  end  of  the  blade.  In  this  manner  the  shear,  moment,  slope, 
and  deflection  at  the  blade  root  may  be  expressed  in  terms  of  the  corresponding  unknown 
tip  values.  The  form  of  Equation  29  that  expresses  the  blade  root  values  in  terms  of  the 


corresponding  unknown  tip  values  is  given  as: 
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I^o]  =  [/^o]  {C?')  +  {®;o) 


(38) 


In  Equation  38,  the  variables  are  as  defined  as  given  in  Equations  30  -  32,  except 

that  the  subscript  0  denotes  the  root  variable  and  the  subscript  T,  the  tip  variable.  The 

column  matrices  of  shear,  moment,  slope,  and  deflection  are  complex,  as  is  the  force  or 
(ato)  column  matrix.  Constants  in  the  4x4  blade  matrix,  [DJ,  are  also  complex.  These 
values  result  from  the  integration  along  the  entire  blade  and  reflect  frequency  as  well  as 

blade  mass  and  stiffness  properties. 

In  the  next  step  we  solve  for  the  unknown  blade  boundary  conditions.  Since  we 

are  dealing  only  with  uncoupled  flatwise  vibrations,  two  boundary  conditions  must  be 

satisfied  at  each  end  of  the  blade.  Application  of  the  known  tip  boundary  conditions 

eliminates  two  unknowns  in  Equation  38  resulting  in  a  2  x  2  [D„  ]  blade  matrix  equation, 
which  is  then  solved  for  the  unknown  tip  values:  (Zj,  67  ). 

For  conventional  free-ended  helicopter  blades,  tip  shears  and  tip  moments  are 

zero.  The  root  boundary  conditions  are  determined  by  the  root  construction  of  the  blade, 

either  articulated  or  rigid  (hingeless).  In  either  case  the  root  vertical  displacement  is  zero, 

but  the  articulated  blade  has  zero  moment  whereas  the  rigid  blade  has  zero  slope  at  the 

root.  Solution  of  the  blade's  forced  harmonic  response  is  then  obtained  by  1)  applying 

appropriate  root  boundary  conditions,  2)  inverting  the  reduced  2x2  complex  blade  matrix 

[DJ,  3)  solving  for  the  unknown  flatwise  tip  slope  and  deflection,  and  4)  using  the 

resulting  4x4  [D]„  to  solve  for  the  shear,  moment,  slope,  and  deflection  along  the  length 

of  the  blade.  Blade,  m  uses  the  2  x  2  matrix  inversion  formula  instead  of  relying  on  the 
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MATLAB’'  function.  This  is  to  reduce  the  possibility  of  numerical  errors  Finally,  this 

procedure  is  repeated  (for  a  total  of  eleven  times)  for  the  steady  (zeroth)  through  tenth 

harmonic  blade  response  with  excitation  applied  successively  at  corresponding  excitation 
frequencies  (0,  Q.  2Q.  3Q .  lOQ).  The  total  response  of  the  blade  at  each  azimuth 

position  is  then  obtained  by  superposition  of  these  results 

C.  OIJTPUT.M 

This  M-file  computes  the  rotor  blade  total  response  and  displays  the  results 
along  with  the  stead\’  and  first  ten  harmonic  responses.  Oiitpiil.ni  also  displavs  graphically 
the  blade  material  properties  These  display  options  are  available  from  a  menu  drix’en 
format  and  the  graphs  may  be  printed  from  the  File  Menu  of  the  figure  Window’". 

The  rotor  blade  total  response  at  each  blade  station  is  calculated  from  the  steady 
and  harmonic  responses  and  is  given  by: 
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(39) 


The  total  response  may  be  viewed  as  a  mesh  plot  for  all  azimuth  angle  positions  or  at 
individual  angles.  The  weight  and  flatwise  stiffness  distribution  as  a  function  of  blade 
radius,  r,  may  also  be  displayed.  Examples  of  sample  output  may  be  viewed  in  the  Results 
and  Recommendations  section  of  this  thesis. 
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IV.  RESULTS 


A.  USER  INSTRUCTIONS 

The  Rotor  Dynamic  Analysis  program  is  a  major  subroutine  of  the  main  JANRAD 
program,  and  can  be  run  after  initiating  the  main  program.  JANRAD  main  program 
installation,  execution,  and  input  requirements  may  be  found  in  Refs,  1  and  2,  Figure  3 
shows  the  main  Execution  Menu  from  where  the  Rotor  Dynamic  Analysis  program  may  be 
accessed.  This  menu  is  accessed  after  either  entering  a  new  data  file  or  editing  a 

***  EXECUTION  MENU  *** 

1 .  Rotor  Performance  Analysis 

2.  Stability  and  Control  Analysis 

3.  Rotor  Dynamic  Analysis 

(Rotor  Performance  Analysis  should  be  run  first) 

4.  Change  data 

5.  Quit 

Enter  a  1, 2,  3,  4,  or  5  » 

Figure  3.  JANRAD  Main  Execution  Menu 

previously  stored  data  file.  It  should  be  noted  that  since  the  Rotor  Dynamic  Analysis 
program  uses  output  data  from  the  Rotor  Performance  Analysis  portion  of  JANRAD,  the 
performance  analysis  should  be  run  first.  If  the  Rotor  Performance  Analysis  has  already 
been  run  for  a  previously  entered  data  file  (and  all  input  parameters  remain  the  same),  and 
a  ".mat"  and  "_p.mat"  already  exist  for  that  file,  the  user  may  proceed  directly  to  the  Rotor 
Dynamic  Analysis  program. 
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The  data  entr>'  and  save  procedures  for  the  Rotor  Dynamic  Analysis  follow  similar 
formats  to  those  of  Refs.  1  and  2.  The  first  screen  of  the  Blade  Dynamic  subroutine 
prompts  the  user  to  either  edit  a  previously  entered  data  file  or  create  a  new  file.  If  the 
edit  file  option  is  chosen,  you  will  be  prompted  to  enter  the  file  name  without  the 
extension.  Just  as  in  the  other  JANRAD  subroutines,  if  the  file  name  is  not  found  the 
Rotor  Dynamic  program  will  allow  the  user  to  try  again  Once  the  file  has  been  loaded, 
the  Blade  Dynamics  edit  menu,  as  shown  in  Figure  4,  will  be  displayed. 

***BLADE  DYNAMICS  EDIT  MENU  *** 

1.  root  boundary  condition 

2.  blade  material  properties 

3.  weight  distribution 

*"  ALL  OTHER  BLADE  INFORMATION  IS  ENTERED  IN  MAIN  JANRAD  MENU  *** 

0.  NO  CHANGES 

Input  the  parameter  to  change; 

Fi}^irc  4.  Blade  Dynamic.s  Edit  Menu 

Whether  the  rotor  blade  input  parameters  are  being  edited  or  a  new  file  created  the 
data  is  entered  in  the  same  manner.  In  the  file  edit  mode,  if  you  choose  not  to  change  the 
value,  pressing  <Enter>  will  re-enter  the  previous  value.  In  the  file  create  mode,  the  user 
will  continue  to  be  prompted  to  enter  a  value,  if  none  has  been  entered,  or  if  a  variable  of 
incorrect  dimensions  is  entered.  For  the  root  boundary  condition,  the  screen  will  prompt 
the  user  to  choose  the  root  boundary  condition  that  corresponds  to  the  construction  of  the 
blade.  The  user  may  enter  the  boundary  conditions  for  an  articulated  rotor  blade.  When 
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entering  the  rotor  blade  material  properties  the  user  has  the  option  of  either  entering  the 
rotor  blade  variable  stiffness  (£/,  Ibsdn"),  or  variable  modulus  of  elasticity  (£,  Ibs/in")  and 
variable  moment  of  inertia  (7^,,  in'*)  separately.  In  either  case  the  screen  will  prompt  the 
user  to  enter  the  row  vector  variable  starting  from  the  tip  of  the  blade  and  ending  with  the 
root;  it  will  also  give  an  example  along  with  the  correct  units.  It  is  noted  that  while  7^. 
typically  varies  with  radius,  in  most  cases.  E  is  constant.  The  program  also  reminds  the 
user  of  the  number  of  elements  the  row  vectors  must  have.  The  number  of  blade  elements 
is  a  variable  that  was  entered  for  the  Rotor  Performance  Analysis  portion  of  JANRAD. 
The  weight  distribution  is  entered  in  the  same  manner.  In  the  file  edit  mode  the  user  is 
returned  to  the  edit  menu  after  editing  any  selection  until  no  further  changes  are  required. 

Once  the  data  file  has  been  edited  or  created  the  program  will  prompt  the  user  for  a 
file  name.  If  the  user  is  in  the  edit  mode,  and  wishes  to  keep  the  same  file  name  they  can 
press  <Enter>  and  the  same  file  name  will  be  used.  Upon  selecting  a  blade  file  name  the 
Rotor  Dynamic  Analysis  program  will  save  the  file  with  a  ".mat"  extension  and  calculate 
the  flatwise  forced  response  of  the  entered  blade. 

B.  OUTPUT 

The  sample  output  is  from  a  helicopter  that  is  based  on  the  Sikorsky  CH-53A  which 
has  a  six  bladed,  articulated  rotor  system.  It  is  taken  from  "Parametric  Investigation  of 
the  Aerodynamic  and  Aeroelastic  Characteristics  of  Articulated  and  Rigid  (Hingeless) 
Helicopter  Rotor  Systems"  [Ref  9].  In  the  report  the  helicopter  is  designated  as  "H2"  and 
its  characteristics  are  shown  in  Appendix  D. 
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Figure  5  shows  the  output  menu  that  appears  after  blade. m  has  finislied  the  required 
calculations.  There  are  four  view  options  from  which  the  user  may  choose.  The 
first  option  (1)  shows  the  steady  and  first  ten  harmonic  responses  of  the  rotor  blade.  They 
are  presented  in  eleven  figure  WINDOWS'  ,  the  first  WINDOW  '  being  the  steady 

CHOOSE  WHICH  OUTPUT  OPTION  YOU  WOULD  LIKE 

1 .  View  the  steady  and  first  ten  harmonic  responses 

2.  View  a  mesh  plot  of  the  flatwise  Shear,  Moment,  Slope  and  Deflection 

at  all  azimuth  positions 

3.  View  the  flatwise  Shear,  Moment,  Slope  and  Shear  at  a  specific 

azimuth  position 

4.  View  the  stiffness  (El)  and  weight  distribution 
0.  Exit 

***  FOR  A  PRINTOUT  CHOOSE  THE  "File"  OPTION  IN  THE  DESIRED  GRAPH  WINDOW  *** 
>> 


Fi};ure  5.  Output  Option  Menu 

response,  the  second  WTNDOW''  is  the  first  harmonic  response,  and  so  on.  For  each 
harmonic  blade  response  there  are  four  subplots  which  showr  blade  flapwise  shear, 
moment,  blade  slope,  and  deflection.  Results  for  the  steady  (or  zeroth)  and  first  ten 
harmonic  responses  of  the  rotor  blade  are  shown  in  Figure  6.  If  option  two  (2)  is  selected, 
the  user  will  see  a  mesh  plot  of  the  total  response  of  the  rotor  blade  over  one 
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Figure  6.  Steady  (2^roth)  Response  of  Sample  Articulated  Rotor  Blade  "H2" 


SHEAR  .  2  HARMONIC  RESPONSE  ^^4MOMENT.  2  HARMONIC  RESPONSE 


(Sai-Nl)  SiN3WOW 


(Nl)  iN3W30\ndSia  1H3A 


Figure  6  (Cont.) .  Second  Harmonic  Response  of  Sample  Articulated  Rotor  Blade  "H2" 
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S'l.- 


Figure  6  (Conf.) .  Third  Harmonic  Response  of  Sample  Articulated  Rotor  Blade  "H2" 
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SHEAR.  3  HARMONIC  RESPONSE  ^  .  ^MOMENT.  3  HARMONIC  RESPONSE 


MOMENT,  4  HARMONIC  RESPONSE 


Figure  6  (Cont.) .  Fourth  Harmonic  Response  of  Sample  Articulated  Rotor  Blade  "H2" 
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SHEAR,  6  HARMONIC  RESPONSE  . MOMENT.  6  HARMONIC  RESPONSE 


Figure  6  (Cont.) .  Sixth  Harmonic  Response  of  Sample  Articulated  Rotor  Blade  "H2" 
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Figure  6  (Cent.)  .  Seventh  Harmonic 
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Figure  6  (Cent.) .  Ninth  Harmonic  Response  of  Sample  Articulated  Rotor  Blade  "H2" 
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Figure  6  (Cont.) .  Tenth  Harmonic  Response  of  Sample  Articulated  Rotor  Blade  "H2" 
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full  revolution  for  the  flatwise  shear,  moment,  slope  and  deflection  (shown  in  Figure  7). 
Option  three  (3),  shown  in  Figure  8,  allows  the  user  to  choose  a  specific  angle  to  view  the 
rotor  blade's  total  response.  Option  four  (4),  shown  in  Figure  9,  allows  viewing  of  a 
graphical  representation  of  the  stifl'ness  modulus  and  weight  distribution. 

In  Figure  6  we  see  the  steady  (zeroth)  response  and  the  first  ten  harmonic  responses 
of  the  sample  articulated  blade  We  note  that  the  dotted  lines  represent  the  imaginary’ 
component  of  the  harmonic  responses.  Consistent  with  complex  notation,  the  "imaginar\'" 
component  of  the  response  leads  the  "real"  component  b\  a  90^  phase  at  the  frequenc\' 
designated.  From  the  steady  (zeroth)  response  it  can  be  seen  that  the  displacement 
represents  the  steady  state  "coning  angle"  of  the  articulated  rotor  system  and  since  there  is 
no  phase  dependency  it  contains  no  imaginary'  component.  The  steady  and  harmonic 
responses  also  confirm  that  all  the  appropriate  boundary  conditions  are  obeyed,  namely: 
that  moments  at  the  blade  hinge  and  tip  are  zero,  deflection  at  the  blade  hinge  is  zero,  and 
vertical  shear  at  the  tip  is  zero.  The  characteristic  mode  shapes  of  the  helicopter  blade  can 
be  seen  to  develop  as  the  frequency  harmonic  increases  up  to  the  third  mode  shape,  typical 
of  the  blade's  tenth  harmonic  response. 


38 


Figure  7.  Mesh  Plot  of  Total  Response  of  Articulated  Rotor  Blade  "H2'’  Over  One  Full  Revolution 


VERT  DISPLACEMENT  (IN) 


Figure  8.  Total  Response  of  Articulated  Rotor  Blade  "H2"  at  a  Chosen  Azimuth  Angle  of  90  Deg. 
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TOTAL  RESPONSE  SHEAR  AT  90  DEG  ^  fl^TAL  RESPONSE  MOMENT  AT 


STIFFNESS  DISTRIBUTION 


Figure  9.  Graph  of  the  Stiffness  Modulus  and  Weight  Dist.  for  Articulated  Rotor  Blade  "H2" 
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BLADE  STATION.  IN 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  modified  Myklestad-Prohl  method  for  rotating  beams  has  been  described  and 
used  in  approximating  solutions  to  aerodynamically  forced  vibrations  of  helicopter  rotor 
blades.  This  type  of  dynamic  response  analysis  has  made  it  possible  to  investigate  the 
dynamic  properties  of  a  helicopter  rotor  blade  of  varying  material  and  sectional  properties, 
JANRAD,  as  a  preliminary  design  tool  has  been  expanded  to  better  meet  the  needs  of  a 
helicopter  design  course.  In  this  respect,  Matlab'*'  has  shown  itself  to  be  the  programming 
language  of  choice  due  to  its  availability,  versatility  in  matrix  manipulation,  and  ease  of 
use.  Although  sample  helicopter  data  is  only  available  for  the  folly  coupled  blade,  results 
show  good  agreement  with  results  given  in  Ref.  9.  Individual  mode  shape  responses  also 
show  good  qualitative  agreement  to  harmonic  excitations. 

B.  RECOMMENDATIONS 

In  this  thesis  we  have  looked  at  the  flatwise  uncoupled  response  of  a  helicopter 
rotor  blade  to  the  effects  of  individual  harmonics  of  airloads.  The  analysis  includes  such 
effects  as  aerodynamic  damping,  and  centrifugal  forces.  Articulated  blade  root  end 
conditions  have  been  treated  in  this  thesis.  An  immediate  improvement  to  this  program 
code  would  be  to  add  hingeless  and  teetering  rotor  blade  root  construction  and  also 
include  provisions  for  blade  lead-lag  or  inplane  response  with  lead-lag  dampers  accounted 
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for  in  the  boundar)'  conditions.  This  would  require  adding  chordwise  dynamics  with  twist 
coupling  into  the  Mykiestad  equations. 

The  Rotor  Performance  Analysis  and  Rotor  Dynamic  Analysis  are  both  based  on  a 
more  extensive  aeroelastic  analysis  that  is  described  in  detail  in  Ref  3.  JANRAD  so  far 
has  looked  at  the  main  cause  of  the  helicopter's  vibrations,  i.e.  the  rotor  system.  For  a 
truly  comprehensive  investigation  of  the  effects  of  rotor  response,  one  should  "marrv'"  the 
rotor  blades  to  the  helicopter  fuselage  As  a  follow-on  thesis  project  the  code  may  be 
expanded  to  include  rotor  blade-fuselage  coupling  through  impedance  matching,  this  is 
discussed  in  detail  in  Ref  3.  Flaving  done  this,  the  forced  response  of  the  helicopter 
fuselage  may  be  calculated  and  we  would  have  available  an  excellent  tool  to  be  applied  in 
both  preliminary  and  detailed  design  of  a  helicopter.  It  could  also  be  used  for  research 
studies  of  rotor/airframe  systems. 
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APPENDIX  A:  DYNAM.M 


%  DYNAM.M 

%  JANRAD:  NPS  Helicopter  Preliminary  Design  Program 
%  Rotor  Blade  Dynamic  Response  Routines 
%  Written  by  Lt  Juan  D.  Cuesta 
%  September  1994 

%  This  program  was  designed  as  an  interactive  preliminary  design  tool 
%  for  rotor  blade  dynamic  analysis  and  design  of  either  an  articulated 
%  or  hingeless  rotor  blade  system.  The  program  provides  the  shear, 

%  moment,  slope  angle,  and  deflection  of  the  flatwise  response  at  any 
%  point  along  the  length  of  a  rotor  blade  to  the  steady  and  first  ten 
%  harmonic  aerodynamic  loads.  This  data  can  then  plotted  at  various 
%  azimuth  blade  positions.  The  input  of  variables  follows  the  same 
%  format  as  written  by  Majs  Bob  Nicholson,  Jr.  and  Walter  Wirth,  Jr. 

%  for  JANRAD. 

%  Variable  List  for  Dynam.m,  Blade. m,  output. m 

%  a  lift  curve  slope 

%  alphaFn  elemental  force  column  matrix 
%  az  azimuth  position  angle 

%  be  root  boundary  condition 

%  cblade2  blade  chord  at  radial  segment,  from  tip  to  root 
%  Cn  flatwise  aerodynamic  damping  on  blade  element 
%  deir  blade  radial  segment  width,  starting  from  blade  hinge 
%  dfn  imaginary  component  of  harmonic  thrust  terms 
%  dFn  real  component  of  harmonic  thrust  terms 
%  dFo  steady  state  thrust  terms 

%  dT  JANRAD  Performance  routine  thrust  output 
%  El  elemental  bending  modulus 

%  En  elemental  modulus  of  elasticity 

%  filenamel  .mat  file  with  janrad  data 
%  filenames  .mat  file  which  contains  blade  data 
%  Fn  distribution  of  thrust  airloads  from  tip  to  root 

%  Ibn  distribution  of  moment  of  inertia  of  blade  elements 

%  Isn  length  of  blade  segment 

%  mn  distribution  of  blade  mass 

%  Mn  flatwise  moment  for  blade  element 

%  naz  number  of  azimuth  sectors 

%  nbe  number  of  blade  elements 

%  omega  rotor  rotational  velocity 

%  omegae  excitation  frequency 
%  Pn  running  transfer  matrix  along  blade  length 
%  psi  azimuth  angle 

%  rho  ambient  air  density 

%  rn  radius,  rotor  blade  radial  segment,  from  blade  hinge 
%  Sn  flatwise  shear  for  blade  element 

%  Thetan  flatwise  slope  of  blade  element 
%  tip_art_bc  articulated  blade  tip  slope  and  deflection  bound,  cond. 

%  tip_rig_bd  hingeless  blade  tip  slope  and  deflection  bound,  cond. 

%  Tn  elemental  radial  tension 

%  Un  Transfer  matrix  between  adjacent  blade  elements 
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%  view  option  variable  for  viewing  choice 
%  Wn  distribution  of  blade  weight 
%  X,Y,Yout  output  data  to  generate  graphs 
%  Yn  flatwise  deflection  of  blade  element 
%  Zn  rotor  blade  element  state  vector 
%  Zroot  rotor  blade  root  state  vector 
%  Ztip  rotor  blade  tip  state  vector 

flag=1; 

flag=exist('filename1'); 
if  flag  ==  0, 
dispC  ') 

dispC  ***  You  must  run  Rotor  Performance  Analysis  first  ***') 
dispf  ***  Press  Any  Key  to  Continue  ***') 

dispC  ') 
pause 
else, 

eval  ([’load  ',filename1  '_p']) 
eval  (['load  ’.filenamel]) 
cic 

dispC  ') 
dispC  ') 

dispC  ***  ROTOR  BLADE  DYNAMIC  ANALYSIS  ROUTINE  ***’) 
dispC  ') 
dispC  ') 

dispC  Do  you  want  to  edit  an  existing  blade  file  or  create  a  new  one?') 

answer3=inputC  1.  edit  existing  file  2.  create  new  file  »'); 

if  answer3==1 , 
cIc 

dispC  ') 
dispC ') 

dispC  ***  load  instructions  *** ') 

dispC  ') 

dispC  A.  Input  the  name  of  the  rotor  blade  data  file  to  edit.') 
dispC  B.  The  file  was  saved  in  your  previous  session’) 
dispC  with  a  ".mat"  extension.') 

dispC  C.  Do  not  include  the  extension  or  quotations. ’) 

dispC  ') 

dispC  ex:  bladef) 

flag=0; 
while  flag  <  1 

filename3=input('  Enter  the  name  of  Blade  data  input  file:  ','s'); 
eval(['flag=existC",filename3,'.mat");']) 
if  flag  <  1 , 
dispC  ') 

dispC  The  file  does  not  exist,  try  again  or  <  Ctrl-C  >') 

dispC  to  exit  program.') 

end 
end 

eval(['load  ',filename3]) 
check4=1 ; 
while  check4  >  0, 
cic 
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***BLADE  DYNAMICS  EDIT  MENU  ***') 


dispC  ') 

dispC 

dispC  ’) 

dispC 

dispC 

dispC 

dispC ') 

dispC ') 

dispC 

dispC  ') 

dispC 


1 .  root  boundary  condition') 

2.  blade  material  properties') 

3.  weight  distribution') 

***  ALL  OTHER  BLADE  INFORMATION  IS  ENTERED  IN  MAIN  JANRAD  MENU  ***') 
0.  NO  CHANGES') 


choicel  =inputC  Input  the  parameter  to  change; '); 
if  isempty(choicel), 
choice1=0; 
end 

if  choicel  ==1, 
cic 

dispC  ’) 

dispCRoot  Boundary  Condition') 
be 

tmp=bc; 
flag=1; 
while  flag  >  0 

bc=inputCRoot  Boundary  Condition  1 .  Articulated  2.  Hingeless  » '); 
if  isempty(bc), 
bc=tmp: 
end 

if  bc==1, 
flag=0; 
elseif  bc==2, 
flag=0; 
else 
dispC ') 

dispC  ***  Enter  a  1  or  2  ***') 
end 
end 
end 

if  choicel  ==2, 
cIc 

dispC  ') 

dispel,  variable  elasticity,  E,  and  variable  moment  of  inertia,  Ixx') 
dispC2.  variable  stiffness,  El') 
option=inputCChoose  1.  OR  2.  »  '); 
if  option==1 , 
cic 

E=En/1e6 

tmp=En; 

disp('1)  Enter  as  a  row  vector  starting  from  the') 

dispC  tip  arid  ending  with  the  root;  ex:  "[18  18.1  ....  21]’”) 

dispC ') 

fprintfC2)  YOU  MUST  ENTER%3.0f  ELEMENTS\n',nbe) 
dispC ') 

disp('3)  Enter  the  modulus  of  ELASTICITY,  E,  distribution  (lbs/in^2  x  10''6): ') 
En=inputC  »').*1e6; 
if  isempty(En), 
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En=tmp; 

end 

while  (length(En)~=nbe), 
disp('l)  Enter  as  a  row  vector  starting  from  the  ') 
disp('  tip  and  ending  with  the  root;  ex:  "(18  18.1  ....  21]"') 
dispC  ') 

fprintf('2)  YOU  MUST  ENTER%3.0f  ELEMENTS\n',nbe) 
dispC  ’) 

disp('3)  Enter  the  modulus  of  ELASTICITY,  E,  distribution  (lbs/in^2  x  10^6):  ') 
En=input('  >>').*1e6; 
end 
cic 
I  bn 

tmp=lbn; 

disp(’l)  Enter  as  a  row  vector  starting  from  the  ') 

disp('  tip  and  ending  with  the  root;  ex;  "[3.9  4.09  ....  15.1]"') 

dispC  ') 

fprintf('2)  YOU  MUST  ENTER%3.0f  ELEMENTS\n',nbe) 
dispC  ’) 

dispC3)  Enter  blade  flapping  moment  of  INERTIA  Ixx.  distribution  (inM): ') 
lbn=inputC  >>'); 
if  isempty(lbn), 
lbn=tmp; 
end 

while  (length(lbn)~=nbe), 
dispCI)  Enter  as  a  row  vector  starting  from  the  ') 
dispC  tip  and  ending  with  the  root;  ex:  "(3.9  4.09  ....  15.1]'") 
dispC ') 

fprintfC2)  YOU  MUST  ENTER%3.0f  ELEMENTS\n',nbe) 
dispC  ') 

dispC3)  Enter  blade  flapping  moment  of  INERTIA,  Ixx,  distribution  (in^4): ') 
lbn=inputC  »'); 
end 

if  existCEr)==1 
clear  El 
end 
end 

if  option==2, 
cIc 
El 

tmp=EI 

dispCI)  Enter  as  a  row  vector  starting  from  the  ') 

dispC  tip  and  ending  with  the  root;  ex:  "[70.2  73,62  ....  271 .8]'") 

dispC ') 

fprintfC2)  YOU  MUST  ENTER%3.0f  ELEMENTS\n',nbe) 
dispC  ') 

dispC3)  Enter  the  bending/stiffness  modulus,  El,  distribution  (lbs  in''2  x  lO'^O):  ') 
EI=inputC  >>').’‘1e6; 
if  isempty(EI), 

El=tmp; 

end 

while  (length(EI)~=nbe), 

dispCI)  Enter  as  a  row  vector  starting  from  the  ') 

dispC  tip  and  ending  with  the  root;  ex:  "[70,2  73.62  ....  271 .8]'") 

dispC  ') 

fprintfC2)  YOU  MUST  ENTER%3.0f  ELEMENTS\n',nbe) 


48 


dispC  ') 

dispCS)  Enter  the  bending  stiffness,  El,  distribution  (ibs  in''2  x  10'^6):  ') 

El=input('  »').*1e6; 
end 

if  exist('En') 
clear  En  Ibn 
end 
end 
end 

if  choice1==3, 
cic 
Wn 

tmp=Wn; 
dispC  ') 

fprintf('l)  YOU  MUST  ENTER%3.0f  ELEMENTS\n',nbe) 
dispC  ') 

disp('2)  Enter  as  a  row  vector  starting  from  the  ’) 

dispC  tip  and  ending  with  the  root;  ex:  "[9.86  9.95  ....  11.96]"') 

dispC ') 

dispC3)  THE  TOTAL  WEIGHT  MUST  BE  GREATER  THAN  THE  AERODYNAMIC) 
fprintfC  PORTION  OF  THE  BLADE:  %6.2f\n’,wblade) 
dispC  ') 

dispCweight  distribution  (Ibs/segment): ') 

Wn=inputC  »'); 
if  isempty(Wn), 

Wn=tmp; 

end 

while  (length(Wn)~=nbe), 
dispC ') 

fprintfCI)  YOU  MUST  ENTER%3.0f  ELEMENTS\n',nbe) 
dispC ') 

dispC2)  Enter  as  a  row  vector  starting  from  the  ') 

dispC  tip  snd  ending  with  the  root;  ex:  "[9.86  9.95  ....  1 1 .96]"') 

dispC ') 

dispC3)  THE  TOTAL  WEIGHT  MUST  BE  GREATER  THAN  THE  AERODYNAMIC) 
fprintfC  PORTION  OF  THE  BLADE:  %6.2f\n',wblade) 
dispC  ') 

dispCweight  distribution  (Ibs/segment): ') 

Wn=inputC  »'); 
end 
end 

if  choice1==0, 
cIc 
dispC ') 
dispC ') 
dispC 
dispC ') 
dispC 
dispC 
dispC 
dispC 
dispC ') 
dispC 
dispC ') 
dispC 
dispC 


***  SAVE  INSTRUCTIONS  ***’) 

A.  Save  the  new  data  to  a  specified  file  name.') 

B.  Do  not  use  an  extension  or  quotations.') 

C.  Use  letter/number  combinations  of  6  characters  or  less.') 

D.  The  file  will  be  saved  with  a  ".mat"  extension.’) 

ex:  bladel') 

E.  If  you  made  no  changes,  press  <  Enter  >,  the  file  will') 
be  saved  with  the  original  name.') 
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flag=1; 
while  flag  >  0 
filename0=filename3; 
filename3=input(’  save  file  as:  '.'s'); 
if  isempty(filename3), 
filename3=filename0; 
end 

clear  filenameO 
if  length(filenamel)  >  6, 
dispC ') 

dispC  use  6  characters  or  less') 
flag=1; 
else 
flag=0; 
end 

eval(['save  ',filename3]) 
check4=0 
end 
end 
end 
end 

%  Creating  a  new  file 
if  answer3==2, 
change=1 ; 
while  change  >  0, 
cic 

bc=input('Root  Boundary  Condition  1.  Articulated  2.  Hingeless  >> '); 
while  isempty(bc), 
dispC ') 

dispCYou  must  enter  a  numerical  value') 

bc=input('Root  Boundary  Condition  1.  Articulated  2.  Hingeless  » '); 
end 
dispC ') 

dispC  Do  you  want  to  enter:') 
dispC ') 

dispC  1-  elasticity,  E,  AND  flapping  moment  distribution,  Ixx') 
dispC  OR') 

dispC  2.  just  the  bending  stiffness,  El,  distribution.') 
clear  optiont 

option1=inputCEnter  1  or  2  >> '); 
while  isempty(optionl), 
optiont  =inputCEnter  1  or  2  »  '); 
end 
cIc 

if  option1==1 , 
dispC ') 

dispCI)  Enter  as  a  row  vector  starting  from  the  ') 
dispC  tip  and  ending  with  the  root;  ex:  "(18  18.1  ....  21]'") 
dispC ') 

fprintfC2)  YOU  MUST  ENTER%3.0f  ELEMENTS\n',nbe) 
dispC  ') 

dispC3)  Enter  the  modulus  of  ELASTICITY,  E,  distribution  (lbs/in''2  x  10''6):') 
En=inputC>>  ').*1e6; 
while  (length(En)~=nbe), 
dispC  ') 

fprintfCI)  YOU  MUST  ENTER%3.0f  ELEMENTS\n',nbe) 
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dispC  ') 

disp('2)  Enter  the  modulus  of  ELASTICITY,  E,  distribution  (lbs/in''2  x  10''6):') 
En=input('»  ').*1e6; 
end 
dispC ') 

disp('1)  Enter  as  a  row  vector  starting  from  the  ') 

disp('  tip  and  ending  with  the  root;  ex:  "[3.9  4.09  ....  15.1]"') 

dispC ') 

fprintf('2)  YOU  MUST  ENTER%3.0f  ELEMENTS\n'.nbe) 
dispC  ') 

dispC3)  Enter  blade  flapping  moment  of  INERTIA,  Ixx,  distribution  (inM): ') 
ibn=inputC»  '); 
while  (length(lbn)~=nbe), 
dispC ') 

fprintfCI)  YOU  MUST  ENTER%3.0f  ELEMENTS\n',nbe) 
dispC  ') 

dispC2)  Enter  blade  flapping  moment  of  INERTIA,  Ixx,  distribution  (inM):  ') 
lbn=inputC»  '); 
end 
end 

if  option1==2, 
dispC ') 

dispel)  Enter  as  a  row  vector  starting  from  the  ') 

dispC  tip  snd  ending  with  the  root;  ex:  "[70.2  73.62  ....  271.8]"') 

dispC  ’) 

fprintfC2)  YOU  MUST  ENTER%3.0f  ELEMENTS\n',nbe) 
dispC ') 

dispC3)  Enter  the  bending  stiffness,  El,  distribution  (lbs  in''2  x  10^6): ') 

EI=inputC»  ').*1e6; 
while  (length(EI)~=nbe), 
dispC ') 

fprintfCI )  YOU  MUST  ENTER%3.0f  ELEMENTS\n',nbe) 
disp('2)  Enter  the  bending  stiffness,  El,  distribution  (lbs  in''2  x  lO'^B): ') 
El=input('»  ').*1e6; 
end 
end 
dispC ') 

disp('l)  Enter  as  a  row  vector  starting  from  the  ') 

dispC  tip  and  ending  with  the  root;  ex:  "[9.86  9.95  ....  1 1 .96]'") 

dispC ') 

fprintf('2)  YOU  MUST  ENTER%3.0f  ELEMENTS\n',nbe) 
dispC  ') 

disp('3)  THE  TOTAL  WEIGHT  MUST  BE  GREATER  THAN  THE  AERODYNAMIC) 
fprintfC  PORTION  OF  THE  BLADE:  %6.2An',wblade) 
dispC  ') 

disp('Enter  weight  distribution  (Ibs/segment)’) 

Wn=input('»  '); 
while  (length(Wn)~=nbe), 
dispC ') 

fprintf('l)  YOU  MUST  ENTER%3.0f  ELEMENTS\n',nbe) 
disp('2)  Enter  weight  distribution  (Ibs/segment)') 

Wn=input('»  '); 
end 
cic 

dispC ') 
dispC ') 
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“  DATA  ENTRY  COMPLETE  ***') 
PLEASE  REVIEW  YOUR  DATA') 


dispC 
dispC 
dispC  ') 

dispC  PRESS  ANY  KEY  TO  CONTINUE') 

pause 

if  option1==1, 

E=En/1e6 

Ibn 

end 

if  option1==2, 

El 

end 

Wn 

dispCDo  you  wish  to  make  any  changes?') 
change=input('0.  No  I.Yes  >>'); 
end 

cic 

dispC ') 
dispC ') 

dispC  ***  SAVE  INSTRUCTIONS  ***') 

dispC  ') 

dispC  A.  Save  the  data  to  a  specified  file  name.') 

dispC  B.  Do  not  use  an  extension  or  quotations. ') 

dispC  S.  Use  letter/number  combinations  of  6  characters  or  less,') 

dispC  D.  The  file  will  be  saved  with  a  ".mat"  extension.') 

dispC ') 

dispC  ex:  bladeT) 
dispC  ') 

dispC  E.  If  you  do  not  enter  a  name,  the  default  is  "bladef" ') 
flag=1; 
while  flag  >  0 

filename3=inputC  save  file  as;  ','s'); 
if  isempty(filename3), 
filename3='blade1': 
end 

if  Iength(filename3)  >  6, 
dispC ') 

dispC  use  6  characters  or  less') 
flag=1; 
else 
flag=0; 
end 
end 

if  existCEI’)==0, 

eval(['save  ',filename3,'  be  En  Ibn  Wn’]) 
end 

if  existCEn')==0, 

eval(['save  ',filename3,'  be  El  Wn']) 
end 
end 
blade 
end 
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APPENDIX  B:  BLADE.M 


%Calculate  steady  state  term 
forj=1:nbe-1, 
dFoG)=(sum(Fn(:,j)))/naz; 
end 


%Caiculate  dFn  (sin  term)  harmonic  matrix 
fori=1;10, 
forj=1:nbe-1, 
dFn(i,j)=0; 
for  k=1:naz, 

Fcum=Fn(k,j)*sin(i*psi(k)/57.3); 

dFn(i,j)=dFn(i,j)+Fcum; 

end 

dFn(i,j)=dFn(i,j)/(naz/2); 

end 

end 


%now  calculate  dfn  (cos  term)  harmonic  matrix 
for  i=1:10, 
forj=1:nbe-1, 
dfn(i,j)=0; 
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for  k=1  :naz, 

fcum=Fn(kJ)*cos(i*psi(k)/57.3); 

dfn(i,j)=dfn(i,j)+fcum; 

end 

dfn(iJ)=dfn(iJ)/(naz/2); 

end 

end 

%  to  iterate  we  switch  the  n's  in  the  myklestad  eqns  (ref  pi  5, 

%  "an  intro  to  helo  dynamics",  E.  R.  Wood)  to  k+l’s 
%  in  the  vector  indices  and  vice  versa 
%  iterations  will  always  be  from  tip  to  root 

dispC ') 

dispC  ***  CALCULATING  THE  AERODYNAMIC  DAMPING  ***') 

%  Calculating  the  flatwise  aerodynamic  damping,  Cn,  on  each  blade  element. 

cblade2=fliplr(cblade): 

for  k=1  :nbe-1 , 

Cn(k)=a*cblade2(k+1)*(rn(k)-rn(k+1))*.5*rho/(12)M*omega*rn(k); 

end 

%Cn=zeros(1 ,23); 

%  the  states  are:  Z=[Sn;  Mn;  Thetan;  Yn] 

%  FIRST  CALCULATE  P  of  equation: 

%  [Z](at  root)  =  [P]*[Z](at  tip)  +  alphaF(at  root) 

dispC  ') 

dispC***  CALCULATING  THE  TRANSFER  MATRICES  FOR  THE  STEADY  AND  10  HARMONIC  ***') 
dispC  ***  BLADE  RESPONSES  ALONG  THE  LENGTH  OF  THE  BLADE  ***’) 

%  CALCULATING  TENSION  first,  since  it  is  independent 
Tn(1)=mn(1)*omega''2*rn(1); 

for  k=1:nbe-1, 

Tn(k+1)=Tn(k)  +  mn(k+1)*omega''2*rn(k+1); 
end 

for  i=0:10; 

omegae=omega*i; 

Pn=eye(4); 
for  k=1  :nbe-1 , 
lsn=rn(k)-rn(k+1); 
if  (existCEI')==0), 

EI=En.*lbn; 

end 

Un=[  1+((mn(k+1)*omegae''2-sqrt(-1)*Cn(k)*omegae)*lsn''3)/(6*EI(k+1)),... 
((mn(k+1)*omegae''2-sqrt(-1)*Cn(k)*omegae)*lsn''2)/(2*EI(k+1)),... 
-(mn(k+1)*omegae''2-sqrt(-1)*Cn(k)*omegae)*(Tn(k)*lsn''3/(6*EI(k+1))+lsn),... 
(mn(k+1)*omegae''2-sqrt(-1)*Cn(k)*omegae); 

lsn+ (Tn  (k)*lsn''3)/(6*E  I  (k+ 1 )) ,  1  +(Tn  (k)*lsn''2)/(2*EI(k+ 1 )) , . . . 

-Tn(k)*(Tn(k)*lsn''3/(6*EI(k+1))+lsn),  0  ; 
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-(lsn''2/(2*EI(k+1))),  -(lsn/EI(k+1)),  (1+Tn(k)*lsn''2/(2*EI{k+1))),  0  ; 
(lsn''3/(6*Ei(k+1))).  (lsn^2/(2*EI(k+1))), -(Tn(k)*lsn''3/(6*EI(k+1))+lsn),  1  ]; 


Unew=Un; 
if  k==1 , 
Ucum=Unew; 
else 

U=[Ucum  Unew]; 
Ucum=U; 
end 


Pn=Un*Pn; 


%  tip  BC  for  all  cases; 

Sn(1)=0; 

Mn(1)=0; 

if  i==0, 

Sn(k+1)=Sn(k)+dFo(k); 

else, 

Sn(k+1)=Sn(k)+dFn(i,k)+sqrt(-1)*dfn(i,k); 

end 

Mn(k+1)=Mn(k)+Sn(k)*lsn; 

Thetan(k+1)=-lsn*Mn(k)/EI(k+1)-lsn''2*Sn(k)/(2*EI(k+1)): 

Yn(k+1)=-lsn^2*Mn(k)/(2*EI(k+1))-lsn^3/(3*EI(k+1)); 

alphaFn(:,k)=[Sn(k+1);  Mn(k+1);  Thetan(k+1);  Yn(k+1)]: 
end 
P=Pn; 

%  SOLVING  FOR  Thetan,  Yn  (NON-ZERO  BOUNDARY  CONDITIONS)  AT  TIP 

%  BOUNDARY  CONDITIONS  ARTICULATED  BLADE  at  root,  Mroot=0,  Yroot=0 
if  bc==1, 

tip_art_bc=-1/(P(2,3)*P(4,4)-P(2,4)*P(4,3))*... 

[P(4,4)  -P(2,4);  -P(4,3)  P(2,3)]*[Mn(nbe);  Yn(nbe)]; 

Ztip=[0;  0;  tip_art_bc]; 
end 

%  BOUNDARY  CONDITIONS  RIGID  BLADE  at  root,  THETAroot=0,  Yroot=0 
if  bc==2, 

tip_rig_bc=-1/(P(3,3)*P(4,4)-P(3,4)*P(4.3))*... 

[P(4,4)  -P(3,4);  -P(4,3)  P(3,3)nThetan(nbe);  Yn(nbe)]; 

Ztip=[0;  0;  tip_rig_bc]; 
end 

Zroot=P*Ztip+alphaFn(:,nbe-1); 

%  CALCULATING  Yout  WHICH  RECORDS  THE  STATES  AT  nbe  STATIONS,  FROM  TIP  TO  ROOT 
Zn=Ztip; 
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for  k=1  :nbe-1 , 
if  k==1 

Zn=Ucum(:,1:4)*Zn; 

Y=Zn+alphaFn(:,k): 

Temp=Y; 

else 

Zn=Ucum(:,4*(k-1)+1 :4*(k-1)+4)*Zn; 
Temp(;,k)=Zn; 
Y(:,k)=Zn+alphaFn(:,k); 
end 
end 

YoLit=[Zfip  Y]: 

Yout2(4*i+1 :4*i+4,:)=Yout; 
Xout=[rn/12]./R; 
end 
output 
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APPENDIX  C:  OUTPUT.M 


%  Subroutine  to  view  options 
view=1 ; 
while  view  >  0, 
dispC ') 

dispC  ***  BLADE  DYNAMICS  OUTPUT  MENU  ***') 

dispC ') 

dispC  CHOOSE  WHICH  OUTPUT  OPTION  YOU  WOULD  LIKE') 

dispC  ') 

dispC  1 .  View  the  steady  and  first  ten  harmonic  responses’) 
dispC ') 

dispC  2.  View  a  mesh  plot  of  the  flatwise  Shear,  Moment,  Slope  and  Deflection ) 
dispC  azimuth  positions') 

dispC  ') 

dispC  3.  View  the  flatwise  Shear,  Moment,  Slope  and  Shear  at  a  specific) 
dispC  azimuth  position') 

dispC  ') 

dispC  4.  View  the  stiffness  (El)  and  weight  distribution') 
dispC  ') 

dispC  0.  Exit') 
dispC  ') 

dispC  ***  FOR  A  PRINTOUT  CHOOSE  THE  "File"  OPTION  IN  THE  DESIRED  GRAPH  WINDOW  ***') 
view=inputC  » '); 
dispC ') 

%  viewing  the  steady  and  the  first  ten  harmonic  responses 
if  view==1 , 
dispC ') 
for  i=0:10, 
figure(i+1) 
subplot(2,2,1) 

plot(Xout,real(Yout2(4‘i+1,:)),'b-',Xout,imag(Yout2(4*i+1,:)),'r-.');grid 
title(sprintfCSHEAR  ,%3,0f  HARMONIC  RESPONSE',!)) 
xlabeK'r/R') 

ylabelCSHEAR  FORCES  (LBS)') 
subplot(2,2,2) 

plot(Xout,real(Yout2(4*i+2,:)),'b-',Xout,imag(Yout2(4*i+2,;)),'r-.');grid 

xlabel('r/R') 

ylabelCMOMENTS  (IN-LBS)’) 

title(sprintf(' . MOMENT,  %3.0f  HARMONIC  RESPONSE',!)) 

subplot(2,2,3) 

plot(Xout,(real(Yout2(4*i+3,:)))*180/pi,'b-’.Xout,(imag(Yout2(4*i+3,:)))*180/pi,'r-.');grid 

title(sprintf('SLOPE,%3.0f  HARMONIC  RESPONSE',!)) 

xlabel('r/R') 

ylabelCSLOPE  (DEG)') 
subplot(2,2,4) 

plot(Xout,real(Yout2(4*i+4,;)),'b-',Xout,imag(Yout2(4*i+4,:)),'r-.’):grid 
title(sprintf('DEFLECTION,%3.0f  HARMONIC  RESPONSE',!)) 
xlabel('r/R') 

ylabelCVERT  DISPLACEMENT  (IN)') 
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end 

end 


%  viewing  the  mesh  plot  for  the  total  response 

if  view==2, 

figure(l) 

for  j=1  :length(psi), 

Yout3=Yout2(1 :4,:);  %steady  state 
for  i=1 :10,  Yoharmonics 

Ynew=real(Yout2(4*i+1 :4*i+4,:)).*sin(psi(j)/57.3*i)+imag(Yout2(4*i+1 :4*i+4,:)).*cos(psi(j)/57.3*i): 
Yout3=Yout3+Ynew; 
end 

YoutS(G-1)+1,:)=Yout3(1,;); 

YoutM(G-1)+1,:)=Yout3(2,:): 

YoutTh((j-1)+1  ,:)=Yout3(3,:); 

YoutY((j-1)+1,:)=Yout3(4,;); 

end 

subplot(2,2,1) 

mesh(Xout,psi,YoutS) 

title(TOTAL  RESPONSE  SHEAR  PLOT) 

xlabel(’r/R') 

ylabel('psi,  deg') 

zlabel('SHEAR  FORCES  (LBS)') 

subplot(2,2,2) 

mesh(Xout,psi,YoutM) 

titleCTOTAL  RESPONSE  MOMENT  PLOT) 

xlabel('r/R') 

ylabel('psi,  deg') 

zlabelCMOMENTS  (IN-LBS)') 

subplot(2,2,3) 

mesh(Xout,psi,YoutTh) 

titleCTOTAL  RESPONSE  SLOPE  PLOT) 

xlabel('r/R') 

ylabelCpsi,  deg') 

zlabeK'SLOPE  (DEG)') 

subplot(2,2,4) 

mesh(Xout,psi,YoutY) 

titleCTOTAL  RESPONSE  DEFLECTION  PLOT') 

xlabel('r/R') 

ylabelCpsi,  deg') 

zlabeK'VERT  DISPLACEMENT  (IN)') 
end 

%  viewing  the  total  response  at  a  specific  azimuth 
if  view==3 
cic 

flag=1; 

pic=0; 

while  flag  >  0, 
pic=pic+1 ; 
dispC ') 

az=input('Enter  the  azimuth  angle  at  which  you  wish  to  see  the  total  response  (deg):  '); 
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Yout3=Yout2(1:4,:);  %steacly  state 
fori=1:10,  %harmonics 

Ynew=real(Yout2(4*i+1:4*i+4,:)).*sin(az/57.3*i)+imag(Yout2{4*i+1:4*i+4,;)).*cos(az/57.3*i); 

Yout3=Yout3+Ynew; 

end 

figure  (pic) 
subplot(2,2,1) 

plot(Xout,Yout3(1  ,:),'b-');grid 

title(sprintf(TOTAL  RESPONSE  SHEAR  AT%3.0f  DEG'.az)) 
xlabel('r/R') 

ylabelCSHEAR  FORCES  (LBS)') 
subplot(2,2,2) 

plot(Xout,Yout3(2,:),'b-');grid 

title(sprintf('TOTAL  RESPONSE  MOMENT  AT%3.0f  DEG’.az)) 
xlabel('r/R') 

ylabelCMOMENTS  (IN-LBS)’) 
subplot(2,2,3) 

plot(Xout,Yout3(3,:)*1 80/pi, ’b-');grid 

title(sprintf(’TOTAL  RESPONSE  SLOPE  AT%3.0f  DEG’.az)) 
xlabel(’r/R’) 

ylabeK’SLOPE  (DEG)') 
subplot(2,2,4) 

plot(Xout,Yout3(4,:),'b-’);grid 

title(sprintf(’TOTAL  RESPONSE  DEFLECTION  AT%3.0f  DEG’.az)) 
xlabel('r/R’) 

ylabeK’VERT  DISPLACEMENT  (IN)') 
dispC  ’) 

disp(’Do  you  want  to  see  another  azimuth  angle?’) 
flag=lnput('  0)  No  1)Yes  »’); 
end 
end 

if  view==4 
cic 

subplot  (2,1,1) 
plot(rn,EI./1e6);grld 

title([’STIFFNESS  DISTRIBUTION  FOR  filenames]) 

xlabel(’BLADE  STATION,  IN') 

ylabel(’FLATWISE  STIFFNESS,  Elxx,  LB*IN''2x10''6’) 

subplot  (2,1 ,2) 
plot(rn,Wn);grid 

title([’WEIGHT  DISTRIBUTION  FOR  ’.filenames]) 
xIabeK’BLADE  STATION,  IN’) 
ylabel(’WEIGHT,  LB/IN’) 
end 
end 
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APPENDIX  D:  H2  SAMPLE  HELICOPTER 


The  following  is  a  listing  of  the  performance  data  file  which  contains  both  input  and 
calculated  parameters  based  on  the  sample  helicopter  "H2",  taken  from  Ref  5. 

filename:  h2_3,prf 

Forward  velocity  =  150  kts 

Temperature  =  59  degs  F 

Pressure  altitude  =  0  ft 

Gross  weight  =  33000  lbs 
Number  of  blades  =  6 

Rotor  radius  =  36.00  ft 
Blade  mean  chord  =  1.97  ft 
Blade  twist  =  -6.00  degs 
Blade  lift  curve  slope  =  5.73 
Blade  weight  =  246.43  lbs 
Rotational  velocity  =  19.37  rads/sec 
Blade  grip  length  =  9.92  ft 
Hinge  offset  =  2.00  ft 
Equivalent  flat  plate  area  =  44.00  ^^^2 
Vertical  projected  area  =  495.00  ft^2 
Wing  area  =  0.00  fR2 
Wing  span  =  0.00  ft 
Wing  CL  =  0.00 
Wing  CDo  =  0.0000 
Wing  efficiency  factor  =  0.00 
Horizontal  tall  area  =  75.00  ft''2 
Horizontal  tail  span  =  15.00  ft 
Horizontal  tail  CL  =  0.40 
Horizontal  tail  CDo  =  0.01 1 5 
Vertical  tail  area  =  84.00  ft''2 
Vertical  tail  span  =  14.00  ft 
Vertical  tail  CL  =  0.40 
Vertical  tail  CDo  =  0.01 14 
Fuselage  drag  =  3356  lbs 
Rotor  drag  =  482  lbs 

Wing  lift  =  0  lbs 

Wing  drag  =  0  lbs 

Horizontal  tail  lift  =  2288  lbs 
Horizontal  tail  drag  =  187  lbs 

Vertical  tail  side  force  =  2563  lbs 
Vertical  tail  drag  =  248  lbs 

Auxiliary  thrust  =  0  lbs 

Tip  path  angle  =  7.90  degs 
Rotor  coning  angle  =  6.65  degs 
Location  of  mean  thrust  (r/R)  =  0.73 
Collective  pitch  at  .7  r/R  =  9.23  degs 
1st  lat  cyclic  term-AI  (deg)  =  2.87 
1st  long  cyclic  term-BI  (deg)  =  -5.75 
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solidity  =  0.105 

CT/sigrna  =  0.063 

CQ/sigma  =  0.0052 

CH/sigma  =  0.0010 

Tip  mach  of  the  adv.  blade  =  0.845 

Advance  ratio  =  0.360 

Rotor  thrust  required  (TPP)  =  31006  lbs 

Rotor  power  required  =  3237  h.p. 

Rotor  torque  =  91908  ft-lbs 

The  following  is  a  listing  of  the  blade  data  file  which  contains  the  material  properties 
of  the  "H2"  articulated  rotor  system,  taken  from  Ref  5  pg  20.3 


RADII 

(IN.) 

CHORD 

(IN.) 

SEGMENT 

WEIGHT 

(LB.) 

FLATWISE 
Ixx  (IN.") 

432.00 

23.65 

6.97 

1,87 

410.40 

23.65 

19.60 

3.75 

388.80 

23.65 

14.51 

4.20 

367.20 

23.65 

15.77 

4.75 

345.60 

23.65 

15.74 

5.05 

324.00 

23.65 

17.14 

5.42 

302.40 

23.65 

10.29 

5.58 

280.80 

23.65 

17.63 

5,68 

259.20 

23.65 

17,02 

5.80 

237.60 

23.65 

18.04 

5.95 

216.00 

23.65 

17.71 

6.05 

194.40 

23.65 

16.61 

6.18 

172.80 

23.65 

16,42 

6.25 

151.20 

23.65 

16.42 

6.40 

129.60 

23,65 

17.33 

6.45 

108.00 

9.71 

18.46 

8.10 

86.40 

9.71 

38.11 

24.10 

64.80 

9.71 

124.62 

100.00 

43.20 

9.71 

77.12 

100.00 

24.00 

9.71 

84.59 

100.00 
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